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Pauli’s exclusion principle has a strong impact on the properties of most fermionic quantum systems. Remark¬ 
ably, the fermionic exchange symmetry implies further constraints on the one-particle picture. By exploiting 
those generalized Pauli constraints we derive a measure which quantifies the influence of the exchange symme¬ 
try beyond Pauli’s exclusion principle. It is based on a geometric hierarchy induced by the exclusion principle 
constraints. We provide a proof of principle by applying our measure to a simple model. In that way, we conclu¬ 
sively confirm the physical relevance of the generalized Pauli constraints and show that the fermionic exchange 
symmetry can have an influence on the one-particle picture beyond Pauli’s exclusion principle. Our findings 
provide a new perspective on fermionic multipartite correlation since our measure allows one to distinguish 
between static and dynamic correlations. 
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I. INTRODUCTION 

The properties and the behavior of most fermionic quan¬ 
tum systems strongly rely on Pauli’s famous exclusion princi- 
pie m. This principle defines a constraint on the one-particle 
picture: For any W-fermion state \^n) the occupancies of 
one-particle states \tp) are restricted, 0 < < 1. 

Indeed, this constrains the one-particle reduced density matrix 
(1RDM) pi = Ntrw-i [|T r j /v)(^ r Ar|], obtained by tracing out 
N — 1 fermion, according to 

0 < pi < 1. (1) 

Equivalent to this matrix relation, the natural occupation num¬ 
bers A^, i.e. the eigenvalues of the 1RDM, are restricted, 

0 < Xi < 1 . ( 2 ) 

These Pauli constraints (PC) play an important role for var¬ 
ious physical phenomena with remarkable consequences for 
both, the micro and the macro world. On a microscopic scale 
they are the basis of the ‘Aufbau principle’ for atoms and nu¬ 
clei. For macroscopic systems the Pauli exclusion principle is 
responsible for the stability of matter mis. In particular, for 
an ideal gas of fermions it implies an effective pressure, the 
Fermi degeneracy pressure, explaining on a rudimentary level 
the stability of neutron stars. Further emerging phenomena 
are the Fermi hole and the Pauli spin blockade. They originate 
from the fact that electrons with parallel spins are forbidden 
to sit at the same position. 

This universal relevance of Pauli’s exclusion principle is 
obvious for weakly correlated systems: All PC are (approxi¬ 
mately) saturated, i.e. one observes for each occupation num¬ 
ber either A^ « 1 or A* ~ 0. Such (approximate) pinning by 
all PC is the typical behavior within the Landau-Fermi theory. 
Even for strongly correlated systems one observes this quasip- 
inning by PC since at least the largest occupation numbers are 
very close to one and the smallest ones are very close to zero. 
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Despite the success on its own, the PC in the one-particle 
picture are de facto a consequence of the antisymmetry of the 
W-fermion wavefunction Ena. This exchange symmetry is 
a much stronger requirement than Pauli’s exclusion principle 
concerning the one-particle picture, only. It has been crucial, 
e.g., for the heuristic approach to the fractional quantum Hall 
effect. Using a Jastrow type wave function involving a product 
of antisymmetric two-particle functions (zi — Zj) p with p an 
odd, positive integer, Laughlin succeeded to explain a huge 
class of plateaus with fractional Hall conductivity ED. 

Since the exchange symmetry concerns the W-particle pic¬ 
ture an important question arises: Does the fermionic ex¬ 
change symmetry have an influence on one-particle properties 
beyond the Pauli exclusion principle? Until a few years ago it 
has even been unclear how to address this question in a mean¬ 
ingful way for concrete systems. This has changed thanks to 
a recent mathematical breakthrough. Motivated by results for 
a very few special cases Q, it was shown that the fermionic 
exchange symmetry in general implies further restrictions on 
the one-particle picture GHna , completing Pauli’s incomplete 
exclusion principle. Those generalized Pauli constraints will 
serve us to eventually address in the present work the long¬ 
standing question above. The corresponding findings exploit¬ 
ing an elegant quantum information theoretical language will 
provide a new perspective on fermionic multipartite correla¬ 
tion since they allow the distinction between static and dy¬ 
namic correlation. 


II. GENERALIZED PAULI CONSTRAINTS 

The so-called generalized Pauli constraints (GPC) take the 
form of linear inequalities 

d 

Dj (A) = > 0 . (3) 

for the natural occupation numbers (NON) A*. These are 
the (decreasingly ordered) eigenvalues of the 1RDM corre¬ 
sponding to the A'-fermion state | ^n}- Here E Z, 
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j = 1,2,..., v( N ' d } < oc and d is the dimension of the un¬ 
derlying one-particle Hilbert space. Geometrically, for each 
fixed pair (TV, d) the family of GPC, together with the nor¬ 
malization and the ordering constraints Ai > ... > A^ > 0, 
defines a (d— 1)-dimensional poly tope V^ N,d>) of allowed vec¬ 
tors A = (A i)f =1 . 

Each of the Dj( A) determines a (d — 2)-dimensional facet 

jr( N i d ) 0 f p(iV,d) } 

= {A € V (N ’ d) I Dj(X) = 0} . (4) 

The remaining facets of V^ N,d ^ are defined by the satura¬ 
tion of the ‘ordering constraints’, i.e. by A^ = A^+i for 
i = 1 , ..., d — 1 . 

To compare the GPC with the PC, notice that the PC 0 
alone define a d-dimensional ‘Pauli hypercube’ C ^ = [0, lp. 
Accounting for the ordering and normalization of the NON, A 
is then restricted to the Pauli simplex 

d 

s (AT,d) = {A 11 > Aa > ... > A d > 0 = N} , (5) 

i=1 

which is a (d — 1)-dimensional subset of the hypercube C^ d \ 
It contains the polytope as a proper subset confirming 

the incompleteness of Pauli’s exclusion principle. For an 
illustration see Fig. |T] In the following we will typically skip 
the indices ‘AT’ and ‘d’. 


III. THE PINNING PHENOMENON 

The potential influence of the fermionic exchange symme¬ 
try beyond Pauli’s exclusion principle is strongly linked to the 
boundary of the polytope corresponding to the saturation of 
a GPC. First, it had been argued EU Q3 that the ground 
state minimization process of the energy expectation value 
(\F|iT|4/) for a model Hamiltonian H may get stuck on the 
boundary of the poly tope V, since any further minimization 
would violate some GPC. This pinning effect by GPC is phys¬ 
ically relevant because it potentially restricts the dynamics of 
the corresponding system whose NON A can never leave the 
polytope 113, 14]. This is in close analogy to the implications 
of the PC which forbid fermions at low temperature to de¬ 
cay to lower lying occupied energy levels. Second, pinning as 
an effect in the one-particle picture allows one to reconstruct 
the structure of the corresponding TV-fermion quantum state 
which in addition is significantly simplified. 

Fet us explain the important latter point. For given (AT, d) 
an arbitrary AT-fermion state |4>) can be expressed as a super¬ 
position of (^) Slater determinants |i) = a| x ... aj |0), built 
from its own natural orbitals \i u ), the eigenstates of the cor¬ 
responding 1RDM. Now, let us assume that given A is pinned 
by a GPC Dj , i.e. Dj( A) = 0. Then, the corresponding su¬ 
perposition is reduced to a subset Zd . containing only those 
configurations i for which Dj\i) = 0 fill . Dj follows from 




FIG. 1. Left side: Schematic illustration of the polytope V, which is 
a proper subset of the light-gray ‘Pauli simplex’ E. The Hartree- 
Fock point is shown as red dot. Right side: For some GPC 
(D i but not D2) a small distance of A to the boundary of the sur¬ 
rounding Pauli simplex enforces a small distance Dj( A) of A to the 
polytope facet Fj . 

Dj (Eq. 0) by replacement of A^ by the occupation operator 
hi of the natural orbital | i). Accordingly, the general superpo¬ 
sition of (£) Slater determinants is reduced to (see also right 
side of Fig. [3]) 

|^(pm)} _ ^ Ci \i) . (6) 

This is quite similar to pinning by a PC, i.e. A^ = 1 or A^ = 0. 
In that case, only those Slater determinants | i) contribute to 
|^) which all do or all do not contain i, respectively. 

Analytical 114 - F7| and numerical investigations [ 18r23| 
have shown that pinning of the corresponding ground state 
NON A to the polytope’s boundary does not occur, in 
general. Whereas pinning has been found for some model 
systems El mm other systems exhibit quasipinning, 
only [12 El [T8l—20ll. In the latter case the vector of NON 
A is close to but not exactly on the polytope’s boundary and 
various implications of pinning hold at least approximately 
(241 . 


IV. HIERARCHY OF PINNING 

A crucial point is that the intersection of the polytope’s 
boundary dV and the boundary <9E of the Pauli simplex is 
not empty. In the following we will show that the set of inter¬ 
section points exhibits a hierarchical structure which in par¬ 
ticular describes the redundancy of the incomplete PC given 
the more restrictive GPC. As a consequence of the inclusion 
relation V C E, already observed but not further elaborated 
in lfl9l . one finds for A G V that A e 9E => A e dV. In 
particular, this applies to the specific part of the boundary of 
E corresponding to pinning by some PC. However, since V 
is not only defined by GPC but also by ordering constraints 
A^ — A^+i > 0, A E dV does not always correspond to pin¬ 
ning by GPC. To explore this in more detail, we express the 
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sets of PC in close analogy to Eq 0 in the compact form 

r d 

S r ,sfy = ^^(1 — A*) + A j > 0 (7) 

i =1 j=d+l—s 

with r < AT, s < d — N, r + s > 0, and introduce the 
corresponding simplicial facets of E 

^r,s = {A G E | 5 r?s (A) = 0} 


= E| Al = ... = Ar = 1 . (8) 

X d+l-s = ■ • • = X d = 0 


Clearly, there is an inclusion relation for those facets: Ejv,d- at 
is a 0-dimensional simplex, coinciding with the Hartree-Fock 
point A hf = ( 1 , •. •, 1 , 0 ,...). E^d-iv is contained in the 
next larger facet, EW-m-jv -2 which is a 1 -dimensional sim¬ 
plex (25). This continues up to the largest two facets Ei j0 and 
Eo,i both of dimension d — 2. In general we have 

E r?s C E r / jS / r > r' A s > s'. (9) 

Similar to A G Fj , pinning by some PC, i.e. A G E r?s , 
implies strong structural simplifications for the correspond¬ 
ing TV-fermion state: r fermions are frozen in the first r 
natural orbitals 11 ),..., |r) and the last 5 natural orbitals 
\d — s + 1 ),..., \d) are inactive, i.e. 

|*> =a|...4|^')- ( 10 ) 

The (TV — r)-fermion state \S&') lives in the corresponding 
active space , |\F') G A N ~ r [h[ d ~ r ~ s ^], where H[ d ~ r ~ s ^ is 
spanned by the active orbitals |r + l),...,|d — s). 

From a geometrical viewpoint, for any pair (r, s) pinning 
by the corresponding PC ([7]) implies pinning to the facet Fj of 
a given GPC Dj if and only if 

PnE^cFj. ( 11 ) 

It is worth noticing that V D E r?s coincides with the poly¬ 
tope V' for the setting (TV', d') M (TV — r,d — r — s) (see 
e.g. Ref. dU). In other words, Eq. (JTTJ) states that VA' G V 
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Based on these observations, natural classes C r?s of GPC 
arise, 

C r , s = {Dj I V n E r , s c Fj} . (13) 

The class C r?s contains exactly those GPC that are saturated 
whenever the corresponding PC ([7]) are saturated. Conse¬ 
quently, the classes C r?s describe the redundancy of the in¬ 
complete PC given the more restrictive GPC. Furthermore, the 
hierarchy <0 implies a hierarchy for these classes, 

C r , s C r<r'As<s', (14) 

and in that sense a partial ordering on {(r, 5)}. By the use 
of a linear program, we determine all classes and calculate 
for each GPC Dj all its minimal (r, s ) with respect to this 
partial ordering, i.e. the smallest pairs (r, s) such that Dj 
still belongs to the corresponding class C r?s . The results 
are listed in the ‘Supplemental Material’ for the three cases 
(N, d) = (3,10), (4,10), (5,10). 



FIG. 2. For the Borland-Dennis setting, (TV, d) — (3,6), the reduced 
polytope of possible spectra (A 4 , A 5 , Ae) is shown (the other three 
NON follow from the equations ( fT6] ). It is spanned by the vertices 
vhf = (0, 0,0), ui = (1/2,1/2,1/2), u 2 = (1/2,1/4,1/4) and 
V 3 = (1/2,1/2,0). GPC < fl7] i takes the form A 5 + A 6 — A 4 > 0 and 
the Pauli Simplex E is spanned by the vertices vhf,v 1 , U 3 and w. 

V. EXAMPLE FOR THE HIERARCHY OF PINNING 

We illustrate the hierarchy of simplicial facets of the Pauli 
simplex E and the induced hierarchy of GPC for the particu¬ 
lar example of the so-called Borland-Dennis setting, (TV, d) = 
(3,6). The corresponding polytope is defined by the con¬ 
straints (7) 

Ai > A 2 > A 3 > A 4 > A 5 > A 6 > 0 (15) 

Ai + A 6 = A 2 + A 5 = A 3 + A 4 = 1 (16) 

D(A) = 2 — (Ai + A 2 + A 4 ) > 0 . (17) 

It is quite special that some GPC take the form of equal¬ 

ities. Due to Eq. © we can choose A 4 ,A5,A6 as the 
independent variables and present the corresponding reduced 
polytope of possible vectors (A 4 ,A 5 ,A 6 ) in Fig. [ 2 ] Also 
various facets of the Pauli simplex are shown: The smallest 
one is the 0-simplex E 3?3 given by the Hartree-Fock point 
(A4, A5, Ae) = ( 0 , 0 , 0 ). It is contained in the next larger facet, 
the 1 -simplex E 2?2 defined by Ai = A 2 = 1 , A 5 = A 6 = 0. 
E 2>2 is then contained in the largest facet, Ei ? i, a 2 -simplex 
shown in gray and defined by A 6 = 0 and A 5 < A 4 < 
Notice, that it does not make sense to consider the other 
facets E r?s with r ^ s since we used the explicit equations 
( f]~6] ) to reduce the poly tope to three dimensions. The only 
GPC which takes the form of a proper inequality is given in 
this reduced description by A 5 + A 6 — A 4 > 0. The facet 
F d corresponding to pinning is spanned by the Hartree-Fock 
point and the vertices V 2 , v 3 and is shown in red. We can now 
determine whether pinning by this GPC is induced by pinning 
by some PC. First, we observe E 3?3 C Fd, i.e. pinning by 
all PC implies pinning of GPC ( [T7] ). Next, we consider the 
smaller pair (r, s ) = (2,2). Clearly, whenever A G V lies in 
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^ 2,2 it also lies in the red facet (actually A has to coincide 
with the Hartree-Fock point). Thus, constraint (17] ) belongs 
also to the smaller class £ 2 , 2 - Does it also belong to the 
smallest class Ci,i? Yes, it belongs also to that class since 
A G PflS i ? i geometrically implies that A lies on the line 
between the Hartree-Fock point and the polytope vertex U 3 . 
This line is an edge of the pinning facet F D . Hence, the 
minimal pair (r, s ) for GPC (17) is (r, s) = (1,1). 


VII. OPERATIONAL MEANING OF Q 

In addition to its main purpose of measuring the influence 
of the fermionic exchange symmetry in the one-particle pic¬ 
ture beyond Pauli’s exclusion principle the Q -parameter has 
also an operational meaning. To elaborate on this in more de¬ 
tail we recall the selection rule ([ 6 | of Slater determinants in 
case of pinning. This remarkable structural simplification of 
A/'-fermion quantum states holds also (approximately) in case 
of quasipinning. To be more precise, one finds [24, 26], 


VI. THE MEASURE 


ot-D(\) < 1 - \\Px D nh < «+ D{ A), (22) 


Since pinning by PC implies pinning by specific GPC ac¬ 
cording to Eq. ( pT] ), the same also holds for quasipinning. Due 
to the flat geometry of polytope facets this implication yields 
linear upper bounds on the l 1 -distance of A to the correspond¬ 
ing polytope facets of the form (see also right side of Fig. [TJ 

disti (A, Fj ) < Cj disti (A, Y r?s ), (18) 

for all (r, s ) with Dj G C r?s . Since we are mainly interested 
in the values Dj (A) we use disti (A, Fj ) oc Dj (A) (see App.jX] 
for more details) to obtain bounds of the form 


Dj( A) < Cj disti (A, Y r?s ). (19) 

Due to the hierarchy ( p~4| ) we consider (T9| ) only for the mini¬ 
mal pairs (r, s ). The optimal (smallest) prefactors Cj for vari¬ 
ous GPC and all their minimal (r, s) are determined by a linear 
program and listed in the ‘Supplemental Material’. 

Equipped with bounds ( [T9] ), we can now decide whether 
given quasipinning by a GPC follows from possible approxi¬ 
mate saturation of PC or whether it goes beyond that. When¬ 
ever Dj( A) is significantly smaller than its upper bound (19) 
the quasipining is stronger than one could expect from pos¬ 
sible quasipinning by PC. Such nontriviality of quasipinning 
by GPC which represents the influence of the exchange sym¬ 
metry beyond the exclusion principle is quantified by the ratio 
of both sides in ( [19] ). This naturally motivates the following 
definition for a measure 


Qj( A) = - log 10 [-Dj(A)/Cjdisti(A,£ r , s )] . (20) 


For GPC with more than one minimal (r, s), i.e. more than one 
independent upper bound (19), we divide Dj( A) in Eq. (|20|) 


by the maximum of all upper bounds. By construction, Qj( A) 
is nonnegative. For given Qj, the minimal distance of A to 
the poly tope boundary is 10 ^ times smaller than one could 
expect from an approximate saturation of PC. Therefore, for 
Qj < 1 quasipinning is rather trivial and for Qj > 2 quite 
nontrivial. For practical applications we define the overall Q- 
parameter, 


Q = max(Qj). (21) 


A ‘Mathematica’ -package implementing all Q-parameters 
can be obtained from the authors or may be found in he 
‘Supplemental Material’. The relation between Q-parameters 
of different settings (TV, d) is discussed in App.jB] 


for some <a_,<a + > 0 , where in general Px is defined as 
the opera tor p rojecting on the space spanned by {| i)}iex- 
Estimate \22\ also holds for PC S r , s ( A) > 0, where the 
corresponding set Xs is given by the configurations i = 
(!,..♦, r, i r + 1 ,.. •, Tjv), r < i v < d — s + 1 , containing all 


1,..., r but none of d + 1 — s,..., d. Such estimates of the 
form ( [ 22 ] ) emphasize the striking relevance that one-particle 


information can have for the description of many fermion 
quantum systems. Combining condition ([22]) for PC S r - :S with 
([22]) for a GPC Dj G C r , s yields 


0 _ ■ 10 "^' (X) < 


i-\\Px Dj nh 

i-\\p Is ^nh 


< 0 + ■ 10 “^' (X) . 


(23) 


To explain the significance of estimate ( f23] >, recall that pin- 
ning by GPC and PC is in a one-to-one correspondence to 
the structure ([ 6 ]) and ( [TO] ), respectively, with corresponding X. 
Since pinning by PC S r , s implies pinning by GPC D G C r?s 
we have X$ r s C Xjj ;j . Therefore, a given 27-fermion quantum 

state with NON A can be written as 


l*> = l^s) + I *D\s) + |*fl> 


(24) 


with \*s) = Pis r ,.\*)> l^vs) = Pt d .\t s J*) and 
= P T <- |'h), where I/j denotes the complement of 


loj■ Estimate (23) then relates the i 2 -weights of D \s) 
and (T'ij). This is illustrated in Fig. [5] On the left, only 
quasipinning by PC is considered which implies that most 
of the L 2 -weight of |\h) is covered by the configurations 
i G Xg r s . However, given additional nontrivial quasipinning 
by a GPC Dj G C r?s , |^) can be further specified since 
almost all weight outside Xs r s needs to lie in Xp. (shown in 
the middle). This means that whenever quasipinning by GPC 
exceeds quasipinning by PC, Eq. ( [23] ) allows one to find much 
more accurate approximations for |4>) than by just freezing 
r electrons in orbitals 11 ),..., |r) and omitting 8 virtual 
orbitals. This latter case of ‘large’ Q includes the specific 
case r = s = 0, i.e. no approximate saturation of any PC 
(illustrated on the right side). Note also that depending on the 
corresponding GPC D', the set of possible configurations 
might be even smaller then the set Xs illustrated on the left 
and in the middle of Fig. [3] 
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FIG. 3. Schematic illustration of how the L 2 -weight of |^) is dis¬ 
tributed over various configurations. Higher weight is indicated by 
darker colors. On the left we use only the information of quasip- 
inning by PC. Stronger implications follow if additional nontrivial 
quasipinning by GPC is taken into account (middle). Right: even if 
none of the PC is (approximately) saturated, strong quasipinning by 
GPC leads to structural simplifications (see also text) 


VIII. PROOF OF PRINCIPLE: APPLICATION TO A 
SIMPLE MODEL 


By applying the Q-parameter to a simple model we provide 
a proof of principle. This will not only emphasize the signif¬ 
icance of our developed concepts but also conclusively con¬ 
firm that the fermionic exchange symmetry has a relevance in 
the one-particle picture beyond Pauli’s 90 years old exclusion 
principle. We consider the same simple model as in Ref. (T5t 

Hn = Jl (§ Vi + j-Z itoZi) + y Y & “ f i) 2 ’ 

i= 1 ' ' l<i<j<N 

(25) 

where x E R n , n E N, Q = diag(cc^,..., andp^ = |V. 

This analytically solvable 7V-Harmonium model has quite 
a long tradition in physics. Yet, its suggested potential signif¬ 
icance for concrete systems as, e.g., for quantum dots (27l is 
here of secondary importance. 

As dimensionless coupling strengths for the different di¬ 
mensions i = 1,... ,n we introduce Ki = Due to a 

duality of NON, observed in Ref. I15U28U30) and proven in 
Ref. ED, we restrict to Ki > 0. We consider the exemplary 
case of N = 3 fermions in n = 2 spatial dimensions. The cor¬ 
responding coupling regime can be parameterized by k = Ki 
and u )2 /and the results for the ground state are shown in 
Fig. [4] There, the solid black line denotes the crossing of the 
ground state and the first excited state and dashed lines indi¬ 
cate crossing of NON which can change the quasipinning be¬ 
havior, as well m. From the left side we learn that quasip¬ 
inning by GPC becomes stronger whenever the coupling k 
decreases and the trap becomes more anisotropic. However, 
the Q-parameter presented on the right shows that the phe¬ 
nomenon of quasipinning is much more involved than previ¬ 
ously appreciated. For instance, the strong quasipinning in the 
regime k - UJ 2 / 0 J 1 ~ 10 is completely trivial, i.e. it is an im¬ 
mediate consequence of quasipinning by PC. It is surrounded 
from both sides by regimes of highly nontrivial quasipinning 
by GPC. Similar results are found also for larger N and larger 
dimensions Eana. This confirms conclusively the physi¬ 
cal relevance of the GPC (recall also Sec. mu and shows that 
the fermionic exchange symmetry can have an influence on 
the one-particle picture beyond Pauli’s exclusion principle for 
concrete systems. 



100 1 10 100 1000 100 1 10 100 1000 


^2 A1 <*>2/^1 

FIG. 4. Quasipinning results for 3-Harmonium ground state in 
two dimensions with external trap frequencies uj\ , UJ 2 and coupling 
strength n. Minimal distance of A to poly tope boundary is 

shown on the left and its nontriviality is quantified on the right (see 
text for details). 


IX. SUMMARY AND CONCLUSION 


We have explored whether the fermionic exchange symme¬ 
try — a property on the 7V-particle level — can have an in¬ 
fluence on the one-particle picture beyond Pauli’s exclusion 
principle. By exploiting a geometric picture, we first ex¬ 
plored and described in the form of a pinning-hierarchy the re¬ 
dundancy of Pauli’s incomplete exclusion principle given the 
more restrictive generalized Pauli constraints (GPC). Based 
on this hierarchy, we succeeded in constructing a measure 
(‘Q-parameter’) which allows one to quantify for concrete 
systems the influence of the GPC beyond the Pauli exclusion 
principle constraints (PC). 

From an operational viewpoint, the Q -parameter has an 
even more important additional significance, since it quanti¬ 
fies possible structural simplifications for TV-fermion quantum 
states beyond the active space size. In particular, this provides 
a new perspective on fermionic multipartite correlation since 
it allows the distinction between static and dynamic correla¬ 
tion. To explain this we first need to discuss the meaning of 
the PC S r:S ( A) > 0 (Eq. ([ 7 ])). The value S' r?s (A) represents 
the weight of the corresponding 7V-fermion quantum state |\F) 
outside the active space obtained by freezing r electrons (in 
the first r natural orbitals) and skipping s virtual states (the 
last s natural orbitals). As a consequence, S' r?s (A) quantifies 
(26l the numerical accuracy of the ‘Complete Active Space 
Self-Consistent Field’-method (see, e.g., the textbook [32]) 
using the corresponding active space. Since so-called ‘dy¬ 
namical correlations’ (see, e.g., Ref. (33)) are attributed to a 
large number of Slater determinants contributing with small 
amplitudes to |4>), the values {5V ?S (A)} measure dynamical 
correlations. Indeed, if S r , s ( A) « 0 for some r, s sufficiently 
close to N and d — TV, respectively, almost all weight of |\F) 
has to lie in a small active space. In contrast to the val¬ 
ues {5V ?S (A)}, the parameters Qj describe how the weight 
of the 7V-fermion quantum state distributes over configura¬ 
tions within the different active spaces. This is illustrated in 
Fig. [3] On the left side (corresponding to a small Qj) all con¬ 
figurations may in principle contribute to |\F), but with small 
amplitudes (dynamic correlations). This is quite different to 
the middle, where ‘large’ Qj implies that the weight shrinks 
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to only a few configurations, reflecting the presence of ‘static 
correlations’ (33j. Since the type of correlations present in a 
given system determines whether a multi-reference or single¬ 
reference approach is required, the Qj -parameters may prove 
useful as an additional tool for this important decision. 

By applying the Q-parameter to a simple model we pro¬ 
vided a proof of principle for the relevance of our new con¬ 
cepts. Our results confirm conclusively the physical relevance 
of the GPC (recall also Sec. [Ill]) and show that the fermionic 
exchange symmetry can have an influence on the one-particle 
picture beyond Pauli’s exclusion principle for concrete sys¬ 
tems. Our findings also suggest that the phenomenon of 
quasipinning by GPC is much more involved than previously 
appreciated. Equipped with our measure Q, the ongoing de¬ 
bate on the relevance of GPC may change significantly after 
reconsidering and reanalyzing lITTl IT^ IT5l 1T814221 l24l [34H37ll . 
Indeed, most of those papers have not contrasted their respec¬ 
tive findings with the less restrictive Pauli exclusion princi¬ 
ple constraints and have therefore confirmed the physical rel¬ 
evance of the generalized Pauli constraints for the wrong rea¬ 
sons. 

Last but not least, our work can be regarded as the starting 
point of a new research avenue based on Klyachko’s break¬ 
through result solving the one-body pure TV-representability 
problem: The systematic analysis and quantification of the in¬ 
fluence of the fermionic exchange symmetry beyond the well- 
established influence of Pauli’s exclusion principle, possibly 
also in the more- particle picture. 
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Appendix A: Distance of A to poly tope facets 

Let us first consider the restriction of the Euclidean space 
R d to the hyperplane V, 

V = {xeR d \x 1 + ...+x d = N}, (Al) 

where N is the particle number of the setting of interest. Be¬ 
low, we will restrict these considerations to proper spectra 
x = {xi)i = i, i.e. to nonnegative and decreasingly-ordered en¬ 
tries Xi. Consider now a hyperplane in V defined by 

E d = {x G V I D(x) = K (0) + K (1) ^1 + ... + K {d) x d = 0} . 

(A2) 


We are interested in deriving an explicit expression for the l 1 - 
distance of y E V to the hyperplane Ed C V. For this, we 
express the distance of y to a set A c V as 

disti (y, A) = sup ({r | B[ r \y) DA = 0}) , (A3) 


where 


B[ r \y) = {xeV\ \\x — y\\i < r} 


(A4) 


is the / 1 -Ball with radius r > 0 around y and ||x||i = 

J2i =i \ x i\ denotes the / 1 -norm. Obviously, B[ r \y) has the 
shape of a cube intersected with the hyperplane x\+.. .+Xd = 
N, i.e. it is in particular a polytope and has therefore finitely 
many extremal points, called vertices. 


For y 0 A it is instructive to think of ( |A3| ) as the radius 
ro that one obtains by further increasing the originally very 
small radius of B^ (y) until one of its points reaches A. For 
A = Ed we observe that in such a process Ed is first reached 
by a vertex of B[ r \y). This follows from the fact that Ed 
is flat. As an elementary exercise one can determine all the 
vertices of B^ (y), following as 


~ y ±r -{ei - ej), 


(A5) 


where i, j = 1 , 2 , ..., d with i < j and ej is a unit vector in 
j-direction. Hence, to determine ro = disti(y, Ed) we need 
to find roots r > 0 of 


D (42a(y)) = D (y) + a n( K 


— I JC-W _ fid-?)) = 0 


(A6) 


for various i, j = 1, 2 ,..., d with i < j and a = ±. ro will 
then be given by the minimum of all those solutions. For a 
consistency check, notice that r = 0 solves ( |A6| ) for y E Ed 
independent of i, j, a and thus disti (y, Ed) = 0. If y 0 Ed 
we have D (jf) -f 0. For i, j with there is no 

solution of (A6). For i, j with we find a solution 

r > 0 whenever a has the correct sign: 


r = 


2 D{y) 


&) - nO) 


(A7) 


Finally, this leads to 


disti (y, E d ) =r 0 = 


2 D(y) 

^tx,max 




(A8) 


where 


A K,„ 


max {| Ki 

i,j 


?!}• 


(A9) 


In particular, the result ( |A8| ) holds for D a GPC and for all 
y E V and therefore in particular also for all y E V C V. By 
restricting the hyperplane Ed to the poly tope V of possible 
spectra A = (A$)f =1 (i.e. in particular to decreasingly ordered 
and nonnegative A^) we obtain a facet Fd of that polytope. 
The ^-di stance of A E V to Fd is identical to the expres¬ 
sion (A8) only for those A whose minimal distance to Ed is 
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attained in the subset Fd C Ed. However, this is always 
the case for A G V sufficiently close to Fd and Ed, respec¬ 
tively. To summarize, whenever A G V is sufficiently strongly 
quasipinned to a facet Fd its l 1 -distance is related to D(X) by 


disti(A ,F d ) 


2D(X) 


(A10) 


For the case that this distance is not small enough, i.e. in 


particular larger than the distance of A to another facet F' of 
V adjacent to F D the right hand side of ( |A10 ) needs to be 
modified by an additional prefactor larger than one. However, 


it is important to keep in mind that the value D( A) is the 
relevant quantity. Still, for didactic reasons it is instructive 
and convenient to talk about the polytope V formed by all 


possible spectra A and relate D (A) to the distance of A to the 
corresponding poly tope facet F D rather than to E D . 


Appendix B: Relation between Q-parameters of different 
settings 

In this section we briefly recall the required concept of trun¬ 
cation and comment on the relation between Q -parameters of 
different settings. We consider two different settings, (TV, d) 
and (TV', d ') = (TV — A AT, d — Ad). Their polytopes V and 
V' are strongly related according to Ref. [34], 


V' = V fl Fi/\N,Ad—AN • (Bl) 


The concept of a truncated quasipinning analysis follows 
from those observations: The analysis of possible quasipin¬ 
ning of given A G V can be simplified by omitting various 
NON very close to 1 and 0 and exploring quasipinning for the 
remaining NON A' in the truncated polytope V'. The minimal 
distance of A to a GPC-facet of V follows then as the mini¬ 
mal distance of A' to a GPC-facet of V' up to a small error. 
This error is given by a linear form in the neglected NON, 
1 - Ai,..., 1 - A A at, A d+ i_Ad+AAr, • • •, A d - Notice, that the 
concept of truncation allows one to perform a quasipinning 
analysis even if the polytope V is not known yet, provided 
that sufficiently many NON are very close to 1 and 0, respec¬ 
tively. This then leads to a truncated setting (TV', d'), whose 
polytope might already be known. 

The concept of a truncated quasipinning analysis was suc¬ 
cessfully applied to various states of the TV-Harmonium ( [25] ) 
studied in this work. There, the one-particle Hilbert space is 
even infinite-dimensional. However, since almost all NON 
are sufficiently close to 0, we could truncate the analysis to 
the settings (3,10), whose GPC are already known. The cor¬ 
responding truncation error was checked to be always smaller 
than various distances of the truncated A' to the GPC-facets of 
V'. 

The concept of truncation also suggests that the Q- 
parameters of different settings are related. One may expect 
the relation 

Q(1 -,...,1-,A',0+ ...,0 +) = Q'(A') ,VA' eV, (B5) 

S -v-' S -v-' 

AN Ad-AN 


This means that intersecting the polytope V of the larger set¬ 
ting with the hyperplane defined by 

Ai = • • • = Xan = 1 A Ad+i-Ad+A7V = • • ■ = Ad = 0 

(B2) 

yields the poy ltop e V' (which is of course still embedded in 
R d ). Relation (B1) means that for every GPC D' of the smaller 
setting there exists at least one corresponding GPC D of the 
larger setting, an extension of D ', i.e. 


AN 


D(X) — 1 — A i) + D r {Xan+i-> 


d 

+ E 

j=d+l — Ad-\-AN 


KjXj . 


> A d-Ad+AN) 

(B3) 


Eq. (|B3|) means nothing else than 


D'(X') = 77(1^1, A', 0^0 ). 

AN d-Ad-\-AN 


(B4) 


For all remaining GPC D of (TV, d) without a p artner GPC 
in (TV', d') their restriction to the hyperplane (B2) leads to an 
inequality in the remaining NON A' which linearly depends on 
the GPC of (TV', d'). As a consequence, those GPC D can be 
neglected whenever the first AN NON are sufficiently close 
to 1 and the last (Ad — AN) NON are sufficiently close to 0. 


where 1 ~, 0 + means to perform a corresponding limit (avoid¬ 
ing possible expressions of the form 0/0). However, since 
some of the individual prefactors in ( [19] ) may vary slightly for 
larger settings, ( |B5| ) may hold in those cases only approxi¬ 
mately. Yet, this would not have any qualitative influence on 
the conclusion of trivial or nontrivial quasipinning. To elab¬ 
orate on this, we consider two GPC D' i n (TV ', d') and D in 


(TV, d) which are related according to Eq. ( |B3| ). Relation ( |B3 ) 
then implies 


Vr, 5 . F) G C r s “w* D G C r -\-AN,s+Ad—AN 5 (B6) 


where C' rs are the classes jl3l for the setting (TV',d') and 
C r?s those for (TV, d). From Eq. ( |B6| ) we conclude that a pair 
(r',s') is minimal (given D' still in C' r , s ,) is equivalent to 
(r, s) m (r' + ATV, s' + Ad — AN) minimal (given D still 
in C r s ). This yields related optimal bounds (recall ([18]) and 


-D'(A') < c , disti(A , ,E^ a ) (B7) 

D(A) < cdisti(A, Sr+AiV.s+Ad-Aiv) • (B8) 


By restricting the bound for D(X) in Eq. (B8) to the hyper¬ 
plane ( |B2| ) yields bound ( |B7| ) up to a potentially different pref- 
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actor. Indeed, we observe 


c = mm 


D( A) 


Xev \disti(A, H r+ AAr, s +Ad-AAr)/ 
( D(X) 


< ^ min - 

xevn'ZAN^d-AN ydisti(A, Y r+ AAr,s+Ad-A;v) 

. ( D'(\‘) \ 

= mm -^- 

X'ev \ disti(A',EJ. a ) J 


(B9) 


where we used Eqs. ( |BT> , ( |B3l > and ( |B4| ) in the second last 
line. 

Whenever the prefactors c and d are different the corre¬ 
sponding Q-parameters (21), Q for D and Q' for D f , do not 
obey Eq. < |B5l >. As a consequence, the overall Q -parameters 
pT) may also not fulfill relation (B5). Yet, this would not 
have any qualitative influence on the conclusion of trivial or 
nontrivial quasipinning. This is based on the fact that c and d 
differ only for a few related GPC D' and D and then only by 
a few percent. 


the second last column empty. For some generalized Pauli 
constraints there is more than one minimal ‘pair’ (but never 
more than two). Each such minimal ‘pair’ gives rise to an 
upper bound of the form 


Dj (A) < Cj disti (A, £ r>a ) VA G V , (C2) 

where E r?s is defined in Eq. (8). By using similar techniques 
as for the derivation of the l 1 2 3 -distance of A to the polytope 
facet Fd (c.f App. © one shows 


disti (A, Y r?s ) = 2 max (£<'-«■ E A,). (C3) 

i= 1 j=d+l — s 

The minimal coefficients Cj in Eq. ( |C2| ) are determined by a 
linear program and are listed in the last columns of the three 
tables in the ‘Supplemental Material’. For any generalized 
Pauli constraint A which does not belong to any class C r ,s 
we establish an elementary upper bound of the form 

A (A) <Ci -1, (C4) 


Appendix C: Explicit form of the Q-parameter 

The generalized Pauli constraints are known so far only for 
settings (TV, d) with d < 10. We present for the three largest 
ones, (TV, d) = (3,10), (4,10), (5,10), the details of the Q- 
parameter and list them in the ‘Supplemental Material’. We 
expressed the generalized Pauli constraints in the canonical 
form 

d 

Dj(\) = K ; o) +5^ K f\ t > o, (CD 

i= 1 


and the corresponding minimal q are listed in the tables below 
as well. 

Let us consider an example: The seventh generalized Pauli 
constraint of the setting (3,10) has two independent upper 
bounds, described by the two minimal ‘pairs’ {r, d +1 — s} = 
{7}, {1,8}. The corresponding prefactors follow as 1 and 
3/4. This means that the optimal bounds on D 7 are given 
by 

D 7 ( A) < 1 • disti (A, £ 0j4 ) and 

D r ( A)<|-dist 1 (A,E ll3 ). (C5) 


where j = 1,2,..., < oo. For each generalized Pauli 

constraint Dj we determine by resorting to a linear program 
their extremal classes C r?s (10), i.e. the minimal pairs (r, s ) 
provided Dj G C r?s . These pairs (r, s) are presented in the 
‘Supplemental Material’, in the second last column of the re¬ 
spective tables. There, we decode these pairs in the form 
{r, d + 1 — s} and also skip possible entries r = 0 and 
d+1 — s = d+ 1. Notice, that a single entry {t} then stands for 
{0, t}ift>N and for {t , d + 1} if t < N. Some generalized 
Pauli constraints do not belong to any class C r?s , i.e. the cor¬ 
responding polytope facet does not contain the Hartree-Fock 
point. This is indicate by leaving the corresponding entry in 


According to the definition (14), the corresponding Q 7 - 
parameter follows as, 


_ DA A) _■ 

max (disti (A, E 0)4 ), f disti (A, Ei )3 )) 

(C6) 

In the same way, we find for the third generalized Pauli con¬ 
straint of the same setting 


<27 (A) = log i 0 


<23 (A) 


log 10 



(C7) 
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{r.d + 1 - s} 

c 

47 

7 

-1 

-3 

-1 

-2 

0 

-3 

2 

1 

0 

2 

{5,6} 

19/20 

48 

7 

-3 

-1 

-2 

-3 

2 

-1 

0 

2 

1 

0 

{3,8} 

2 

49 

7 

-1 

-2 

-3 

-1 

0 

-3 

2 

1 

0 

2 

{3,8} 

2 

50 

7 

-1 

-3 

-2 

-3 

2 

-1 

1 

0 

0 

2 

{5,6} 

19/20 

51 

7 

-3 

-1 

-1 

-2 

0 

-3 

2 

1 

2 

0 

{5,6} 

19/20 

52 

7 

-3 

-1 

-2 

-1 

0 

-3 

2 

2 

1 

0 

{3,8} 

2 

53 

7 

-1 

-2 

-3 

-3 

2 

-1 

0 

1 

0 

2 

{3,8} 

2 

54 

7 

-3 

-2 

-3 

-1 

2 

1 

0 

-1 

2 

0 

{3, 10} 

2 

55 

7 

-1 

-3 

0 

-1 

-2 

-3 

0 

2 

1 

2 

{U8} 

2 

56 

7 

-1 

-1 

-2 

-3 

0 

-3 

0 

2 

1 

2 

{1,8} 

2 

57 

7 

-3 

-2 

-3 

-1 

2 

-1 

2 

1 

0 

0 

{3, 10} 

2 

58 

7 

-1 

-1 

-2 

-3 

0 

-3 

1 

0 

2 

2 

{5,6} 

31/30 

59 

7 

-3 

-3 

-1 

-2 

2 

-1 

2 

1 

0 

0 

{5,6} 

31/30 

60 

9 

-3 

-4 

-5 

-1 

4 

3 

2 

1 

0 

-2 

{5,6} 

2 

61 

9 

-3 

-4 

-5 

3 

2 

1 

0 

-1 

4 

-2 


11/3 

62 

9 

1 

-5 

0 

-1 

-2 

-3 

-4 

4 

3 

2 


11/3 

63 

9 

1 

-1 

-2 

-3 

-4 

-5 

0 

4 

3 

2 

{5,6} 

2 

64 

4 

-1 

-2 

0 

-2 

1 

0 

-1 

2 

1 

0 

{3,8} 

7/4 

65 

4 

0 

-1 

-2 

0 

-1 

-2 

1 

0 

2 

1 

{3,8} 

7/4 

66 

4 

-2 

-2 

-1 

1 

1 

0 

0 

-1 

2 

0 


3/2 

67 

4 

0 

-1 

-2 

-2 

1 

0 

0 

-1 

1 

2 

{5,6} 

7/10 

68 

4 

-2 

0 

-1 

-1 

0 

-2 

2 

1 

1 

0 

{5,6} 

13/20 

69 

4 

0 

-1 

-2 

-2 

1 

0 

-1 

0 

2 

1 

{3,8} 

7/4 

70 

4 

0 

-2 

-1 

-2 

1 

0 

0 

-1 

2 

1 

{5,6} 

7/10 

71 

4 

-2 

-1 

0 

0 

-1 

-2 

2 

1 

1 

0 

{5,6} 

11/16 

72 

4 

-2 

-2 

-1 

-1 

2 

1 

1 

0 

0 

0 

{5,6} 

4/5 

73 

4 

0 

-1 

-1 

-2 

0 

-2 

1 

0 

2 

1 

{5,6} 

7/10 

74 

4 

-1 

-2 

0 

-1 

0 

-2 

2 

1 

1 

0 

{5,6} 

13/20 

75 

4 

-1 

-2 

0 

0 

-1 

-2 

1 

2 

1 

0 

{3,8} 

7/4 

76 

4 

0 

-2 

0 

-1 

-1 

-2 

0 

2 

1 

1 

{1,8} 

3/2 

77 

4 

0 

-1 

-1 

-2 

0 

-2 

0 

2 

1 

1 

{1,8} 

3/2 

78 

4 

0 

-2 

0 

-2 

0 

-1 

-1 

2 

1 

1 

{1,8} 

3/2 

79 

2 

0 

-2 

1 

0 

0 

-1 

-1 

1 

2 

2 


3/2 

80 

2 

0 

-1 

-2 

1 

0 

-1 

2 

0 

2 

1 

{3,8} 

7/4 

81 

2 

-1 

-2 

0 

-1 

2 

1 

0 

2 

1 

0 

{3,8} 

7/4 

82 

2 

-2 

-1 

1 

0 

0 

-1 

2 

2 

1 

0 

{5,6} 

7/10 

83 

2 

0 

-1 

-1 

-2 

2 

0 

1 

1 

0 

2 

{5,6} 

13/20 

84 

2 

-1 

-2 

0 

1 

0 

-1 

2 

2 

1 

0 

{3,8} 

7/4 

85 

2 

-1 

-2 

1 

0 

0 

-1 

2 

1 

2 

0 

{5,6} 

7/10 

86 

2 

0 

-1 

-1 

-2 

2 

1 

0 

0 

1 

2 

{5,6} 

11/16 

87 

2 

-1 

-2 

0 

-1 

2 

0 

2 

1 

1 

0 

{5,6} 

7/10 

88 

2 

0 

-1 

-1 

-2 

2 

0 

1 

0 

2 

1 

{5,6} 

13/20 

89 

2 

0 

-1 

-2 

-1 

2 

1 

0 

0 

2 

1 

{3,8} 

7/4 

90 

2 

-1 

-1 

-2 

1 

1 

0 

2 

0 

2 

0 

{3, 10} 

3/2 

91 

2 

0 

0 

0 

-1 

-1 

-2 

1 

1 

2 

2 

{5,6} 

4/5 

92 

2 

-1 

-1 

-2 

0 

2 

0 

2 

1 

1 

0 

{3, 10} 

3/2 


















